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Abstract. We describe various aspects of the Al-Salam-Carlitz g-Charlier polyno- 
mials. These include combinatorial descriptions of the moments, the orthogonality 
relation, and the linearization coefficients. 



1. Introduction. 

The Charlier polynomials C®(x) are well-known analytically [4], and have been 
studied combinatorially by various authors [8], [12], [16], [17], [20]. The moments 
for the measure of these orthogonal polynomials are 



;i.l) n n = Y,S{n,k)a h 



k=i 



where S(n, k) are the Stirling numbers of the second kind. The purpose of this paper 
is to study combinatorially an appropriate (/-analogue of C%(x), whose moments are 
a g-Stirling version of (1.1). While studying these polynomials, we use statistics on 
set partitions which are g-Stirling distributed. 

Our main result (Theorem 3) is the combinatorial proof of the linearization 
coefficients for these polynomials. In the q = 1 case, the linearization coefficients 
are given as a polynomial in a, whose coefficients are quotients of factorials (see 
(4.4)). This has a simple combinatorial explanation. However, in the q-c&se the 
coefficients are not the analogous quotients of ^-factorials. They are alternating 
sums of quotients of g-factorials, and thus a combinatorial explanation is much 
more difficult. From the combinatorial interpretations of the polynomials and their 
moments, in terms of weighted partial permutations and set partitions, we deduce a 
combinatorial interpretation for the linearization coefficients of a product of three q- 
Charlier polynomials. We then apply a weight-preserving sign-reversing involution 
defined in five steps. Theorem 3 is obtained by enumerating the remaining fixed 
points. Some of the steps of the involution are quite straight-forward, but some 
others are more complicated. They use more sophisticated techniques such as 
encoding of permutations or set partitions into 0-1 tableaux (cf [6], [18]), which 
are fillings of Ferrers diagrams with O's and l's such that there is exactly one 1 
in each column. They also use interpolating statistics on set partitions, as were 
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introduced by White in [22]. Indeed, the characterization of the final set of fixed 
points uses a bijection \I/g of White between interpolating statistics, making their 
enumeration all the more complicated. 

It turns out that our g-Charlier polynomials are not what have classically been 
called (/-Charlier; in fact they are rescaled versions of the Al-Salam-Carlitz polyno- 
mials [4, p. 196]. Some comparisons to the classical g-Charlier are given in §7. Zeng 
[24] has also studied both families of polynomials from the associated continued 
fractions. 

The basic combinatorial interpretation of the polynomials is given in Theorem 
1. Several facts about the polynomials can be proven combinatorially. The combi- 
natorics of set partitions, restricted growth functions and 0-1 tableaux is discussed 
in §3, and the statistic for the moments is given in Theorem 2. In §4, we state 
our main theorem, Theorem 3, giving the linearization coefficient for a product of 
three g-Charlier polynomials, and we set up the general combinatorial context for 
its demonstration. The five steps of the weight-preserving sign-reversing involu- 
tion proving Theorem 3 are given in §5, and the combinatorial evaluation of the 
remaining fixed points is the subject of §6. 

We use the standard notation for (/-binomial coefficients and shifted factorials 
found in [11]. We will also need 



2. The g-Charlier polynomials. 

We define the g-Charlier polynomials by the three term recurrence relation 



where C-i(x, a; q) = and Cq(x, a; q) = 1. 

It is not hard to show that these polynomials are rescaled versions of the Al- 
Salam Carlitz polynomials [4, p. 196] 

(2.2) C n (x, a; q) = a n U n (- - -tt^t, ^f^)- 

a a(l — q) a(l — q) 

Since the generating function of the U n (x, b) is known [4], we see that 




and 



[n]\ q = [n] q [n-l] q - ■■[!},. 



(2.1) 



C n+1 (x,a;q) = (x - aq n - [n] q )C n (x, a; q) - a[n] q q n 1 C n - 1 (x,a;q), 



(2.3) 




This gives the explicit formula 
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Clearly, we want a q- version of [16], which gives the Charlier polynomials as 
a generating function of weighted partial permutations, i.e. pairs (S,cr), where 
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B C {1, 2, • • • , n} = [n], and a is a permutation on [n] — B. Thus we need only 
interpret the individual terms in (2.4) for a combinatorial interpretation. The inside 
product can be expanded in terms of the (/-Stirling numbers of the first kind. We 
let cyc(a) be the number of cycles of a permutation a and inv(a) be the number of 
inversions of a written as a product of disjoint cycles (increasing minima, minima 
first in a cycle). 

fc-i 

I]> -[*]<?)= E (_l)k-cyc{v) ql nv{v) x cyc{v)^ 

For the sum over k in (2.4), we sum over all (n — k) subsets B C [n]. Let 

inv(B) = J2(b-l), 
beB 

so that the generating function for these subsets is 

(n-k\ 
A 2 /. 



We have established the following theorem. 
Theorem 1. The q-Charlier polynomials are given by 

C n (x,a;q)= J2 q inv ^ )+inv{B \-l) n - cyc{a) a^x cyc{a \ 

BC[n] ae & n _ B 

BC[n] ae Q n _ B 

A combinatorial proof of the three-term recurrence relation (2.1) can be given 
using Theorem 1. An involution is necessary. For more details, we refer the reader 
to [5]. 

3. The moments. 

An explicit measure for the g-Charlier polynomials is known, [4, p. 196]. It is 
not hard to find the n th moment of this measure explicitly. The result is a perfect 
g-analogue of (1.1) 

n 

(3.1) ^ n = J2S q (n,k)a k , 

k=i 

where S q (n, k) is the g-Stirling number of the second kind, given by the recurrence 

(3.2) S q (n, k) = S q {n - 1, k - 1) + [k] q S q (n - 1, k), 
where S q (0, k) = Sq^- In fact, one sees that [13] 



^ n—k 

(3.3) S,(n,k)= —— — k 
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(k + j) 


~k+j~ 




j 



(-iy. 



Clearly (3.1) suggests that there is some statistic on set partitions, whose gener- 
ating function is \i n . This statistic, rs, arises from the Viennot theory of Motzkin 
paths associated with the three-term recurrence (2.1) [20]. We do not give the 
details of the construction here. 

However, let us review some combinatorial facts about (/-Stirling numbers. Set 
partitions of [n] = {1, 2, . . . , n} can be encoded as restricted growth functions (or 
RG-functions) as follow: if the blocks of n are ordered by increasing minima, the 
RG-function w = w\W2 ■ ■ ■ w n is the word such that Wi is the block where i is 
located. For example, if tt = 147|28|3|569, w = 123144124. Note that set partitions 
on any set A can be encoded as RG-functions as long as A is a totally ordered set. 

In [21], Wachs and White investigated four natural statistics on set partitions, 
called Is, lb, rs and rb. They are defined as follow: 

n 

ls(w) = \ {j : j < Wi,j appears to the left of position i}\, 
i=i 

n 

lb(w) = \ {j : j > Wi,j appears to the left of position i}\, 

n 

rs(w) = \{j : j < Wi,j appears to the right of position i}\, 
i=i 

n 

ls(w) = \{j : j > Wi,j appears to the right of position 

i=i 

Thus in the example, ls(n) = 13, lb(n) = 7, rs(n) = 7 and rb(n) = 11. They 
showed, using combinatorial methods, that each had the same distribution (up to 
a constant) on the set RG(n, k) of all restricted growth functions of length n and 
maximum k, and that their generating function was indeed S q (n, k) for rs and lb 

(respectively q( 2 ' S q (n,k) for Is and rb). 

We also use another encoding of set partitions in terms of 0-1 tableaux. A 0-1 
tableau is a pair <p = (A, /) where A = (Ai > A2 > . . . > Afc) is a partition of 
an integer m = |A| and / = (fij)i<j<\i is a "filling" of the corresponding Ferrers 
diagram of shape A with 0's and l's such that there is exactly one 1 in each column. 
0-1 tableaux were introduced by Leroux in [18] to establish a (/-log concavity result 
conjectured by Butler [3] for Stirling numbers of the second kind. 

There is a natural correspondence between set partitions tt of [n] with k blocks 
and 0-1 tableaux with n — k columns of length less than or equal to k. Simply 
write the RG-function w = W1W2 ■ ■ .w n associated to n as a k x n matrix, with a 
1 in position if Wj = i, and elsewhere. The resulting matrix is row-reduced 
echelon, of rank k, with exactly one 1 in each column. A 0-1 tableau (in the third 
quadrant) is then obtained by removing all the pivot columns and the 0's that lie 
on the left of a 1 on a pivot column. Figure 1 illustrates these manipulations for 
tt = 1247|39(12)|568(11)|(10). 

We define two statistics on 0-1 tableaux cp: first, the inversion number, inv(<p), 
which is equal to the number of 0's below a 1 in <p; and the non-inversion number, 
nin(ip), which is equal to the number of 0's above a 1 in if. For example, for ip 
in Figure 1, inv(ip) = 7 and nin(ip) = 8. Note that an easy involution on the 
columns of 0-1 tableaux sends the inversion number to the non-inversion number 
and vice-versa. We call this map the symmetry involution. 
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Zs(tt) = 
lb(n) = 
rs{ji) = 
rbiiz) = 
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Figure 1: Correspondence between partitions and 0-1 tableaux 



It is not hard to see that the inversion number (respectively non-inversion num- 
ber) on 0-1 tableaux corresponds to the statistic lb (resp. Is — ( 2 )) on set partitions. 

Similarly, permutations a of [n] in k cycles can be encoded as 0-1 tableaux with 
n — k columns of distinct lengths less than or equal to n — 1. The correspondence is 
defined by recurrence on n. Suppose a is written as a standard product of cycles. 
If n = 1, then a = (1) corresponds to the empty 0-1 tableau ip = 0. Otherwise, 
let a G & n +i and let <p denote the 0-1 tableau associated to the permutation a in 
which (n + 1) has been erased. There are two cases. If (n + 1) is the minimum of 
a cycle in a, then a corresponds to <p. If (n + 1) is not the minimum of a cycle, 
then it appears in a at a certain position i, 2 < i < n + 1. The permutation a 
then corresponds to the 0-1 tableau <p> plus a column of length n with a 1 in the 
(i — l)-th position (from top to bottom). For example, a = (1, 3, 4, 7, 2)(5, 6)(8) 
corresponds to the following 0-1 tableau. 
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Figure 2: Correspondence between permutations and 0-1 tableaux 



It is not hard to see that under this transformation, the inversion number on 
0-1 tableaux corresponds to the inversion number on permutations, as defined in 
§2. Thus, their generating functions are the (/-Stirling numbers of the first kind 
c q (n, k). 

In [6] , de Medicis and Leroux investigated q and p, (/-Stirling numbers from the 
point of view of the unified 0-1 tableau approach. In particular, they proved com- 
binatorially or algebraically a number of identities involving (/-Stirling numbers. 

For the combinatorial interpretation of the moments of the (/-Charlier polyno- 
mials in terms of set partitions n, we need two statistics. The number of blocks 
^blocks(ji) is one, and the other statistic is rs(jv). 
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Theorem 2. The n th moment for the q-Charlier polynomials is given by 

lj n = a #Wocfcs(7 V s(7r) . 

-n-eP(n) 

As we mentioned, many other (/-Stirling distributed statistics have been found 
[21]. It is surprising that the Viennot theory naturally gives a so-called "hard" 
statistic (rs), not an easy one (e.g. lb, [21]). Other variations on the rs-statistic 
can be given from the Motzkin paths, although the /^-statistic is not among them. 
It can be derived from the Motzkin paths associated with the "odd" polynomials 
for (2.1). 

4. The orthogonality relation and the linearization of products. 

Let L be the linear functional on polynomials that corresponds to integrating 
with respect to the measure for the Charlier polynomials. The orthogonality rela- 
tion is 

(4.1) L{C a n {x)C a m {x)) = a n nl 5 m , n . 
The q- version of (4.1) is 

(4.2) L q (C n (x, a; q)C m (x, a; q)) = a n q(^ [n]\ q 8 m , n . 

Since the polynomials C n (x, a; q) and L q have combinatorial definitions from The- 
orems 1 and 2, it is possible to restate (4.2) as a combinatorial problem. We will 
give an involution which then proves (4.2) in this framework. 

A more general question is to find L(C^ l (x)C!^ 2 (x) ■ ■ ■ C!^ k (x)) for any k. A 
solution is equivalent to finding the coefficients a Uk in the expansion 

C a ni (x)C a n2 (X) ■ ■ ■ (% k _ x (X) = <*n k C a nk ( X) . 

This had been done bijectively for some classes of Sheffer orthogonal polynomials in 
[5], [7], [9], [10]. Moreover, in the g-case of Hermite polynomials, some remarkable 
consequences have been found [15]. 

For the Charlier polynomials, it is easy to see that 
(4.3) 

E HC^(x)C^(x) ■ ■■C a nk {x))±- ■ ■ ■ -a- = e -(e a (*i.-.**)+-+«*(*i.-,t fc )) j 

ni,--,n fe =0 

where is the elementary symmetric function of degree i, [19]. In this case 
L(C^ l C^ 2 • • 'C% ) is a polynomial in a with positive integer coefficients; a combi- 
natorial interpretation of this coefficient has been given ([12] and [23]). For k = 3, 

(4.3) is equivalent to 
(4.4) 

L(ni+n 2 -n 3 )/2j n3+z ( ( ( 

L(C a ni (x)C a n2 (x)CZ 3 (x)) = g /!(n 3 - n 2 + l)l(n 3 - nT+l)^ +n 2 -n 3 -2l)V 

One can hope that L q (C ni (x)C n2 (x)C n3 (x)) is simply a weighted version, with 
an appropriate statistic, of the q = 1 case. However this is false. For example, 

L,(C 2 (x)C 2 (x)C7i(x)) = q(q 2 + 2q + l)a 2 + q(q 3 + q 2 - q - l)a 3 . 

Nonetheless, we have an exact formula for L q (C ni (x, a, q)C n2 (x, a, q)C n3 (x, a, q)), 
which is equivalent to one of Al-Salam-Verma [1]. 

6 



Theorem 3. Let n 3 > ri\ > n 2 > 0. Then 

L q (C ni (x)C n2 (x)C n3 (x)) 
(4.5) [n 3 ]l q 



ni+ri2—ns I 
Z=0 j=0 



[n!-l]l q 



n 2 
1-3 



n\ 




n 2 -l + j 


_ 3 _ 


q 


_n 3 -n 1 +j _ 



3% 



q [n 3 - n 2 + l]\ q 



ni + n 2 - n 3 - I 
3 



where 



+ (n 3 - m +j)(n 3 -n 2 + l) + j(n 3 - n 2 +/). 



The generating function of L q (C ni (x)C n2 (x)C n3 (x)) can be evaluated from The- 
orem 3, yielding 



L q {C ni {x)C n2 {x)C nz {x))j± 2 



+"3 

l 3 



ni,n 2 ,n 3 

(4.6) 



[n 2 }l q [n 3 ]! c 



(-^3; 9)00 (-0*1*2(1 - q);q)oo i4>\ 



at\{l - q), at 2 (l - q) 
-at\t 2 {l - q) 



Letting q — > 1 in (4.6) gives back (4.3) for k = 3. This generating function can also 
be evaluated directly using the measure ([4, p. 196]), the generating function (2.3) 
for the polynomials and a 3 <p 2 transformation. 

More generally, for k > 4, the generating function of L q (C ni (x) . . .C nk (x)) can 
be expressed as a difference of two basic hypergeometric series. This has been done 
by Ismail and Stanton [14] for the Al-Salam Carlitz polynomials, so an equivalent 
formula can be deduced for the g-Charlier polynomials using (2.2). 

Let us set up the combinatorial context in which Theorem 3 will be proven. We 
first introduce notations and conventions that will be used throughout the proof. 
Define 

L q (n u n 2l n 3 ) = {((B h cr^); 7r) = ((£1,01), (B 2 ,a 2 ), (B 3 , 03); n) | 

(Bi, (Ji) is a partial permutation on the set {i} x [rii], 
and 7T is a partition on the cycles of a±, a 2 and a 3 }. 



We will say that an element of the set {i} x [rii] is of color i. When giving exam- 
ples of elements of L q (ni,n 2 ,n 3 ), to simplify notation, pairs (l,i), (2, 1) and (3, i) 
will always be denoted i, i and i respectively. Thus a typical element of L q (8, 7, 10) 
would be described in the following way: B\ = {2, 3}, B 2 = 0, B 3 — {5, 9, 10}, and 

= (I, 5, 7) (8) (T) I (4) I (6) (3, 5) (3, 7) | (1, 4, 2) (6, 7) | (2, 8, 4, 6) (the underlying permu- 
tations (Ji = (1,5, 7) (4) (6) (8), a 2 = (1,4,2)(3,5)(6,7) and a 3 = (T) (2, 8, 4, 6) (3, 7) 
can be recovered from n). 

Note that the lexicographic order on pairs induces a total order on the 

cycles of o"i, a 2 and a 3 , according to their minima. Therefore we can talk about RG- 
functions. We will always use the letter w to denote the RG-function associated to 
7r. In the above example, w = 1231434153. The first cyc(ai) letters of w correspond 
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to the positions of cycles of color 1 in tt, the next cyc(a 2 ) to the positions of cycles 
of color 2, and the last cyc(<j 3 ) letters to the positions of cycles of color 3. We will 
denote by w a , Wf, and w c respectively these portions of w. In the above example, 
we have w a = 1231, Wb = 434, w c = 153, and w = w a WbW c , the concatenation of 
words w a , Wb and w c . 

Finally, we will use the notation Supp(w) (or Supp(a) or Supp(TTi)) to denote 
the underlying set of letters of a word w (or a permutation a or a block tt{ of a 
partition tt respectively). 

From Theorems 1 and 2, we deduce that 

(4.7) L q (C ni (x)C n2 (x)C n3 (x)) = J2 w,((Bi,<7i);7r), 

((B i ,cr i );7r)6L (3 (ni,n2,n3) 

where 

(4.8) u q ((Bi, a,); tt) = u q {B u a^Bi, a 2 )u q (B 3 , a 3 )q rs ^a* blocks ^ , 

and u q (B,a) was defined in Theorem 1, clS ct SI gned monomial in the variables a 
and q. This gives a combinatorial interpretation of the right-hand side of (4.5). 

For q = 1, the negative coefficients of a are counterbalanced by the positive 
coefficients of a, and (4.7) is a polynomial with positive coefficients. Indeed, in 
that case, it is not hard to find a weight-preserving sign-reversing involution on 
L q (ni, n 2 , n 3 ) (cf [5]) whose fixed points ((Bi, o~i); tt) are characterized by 

i) Bi = and = Identity, for i = 1, 2, 3; 

ii) the word w a (respectively w/& and iy c ) contains all distinct letters, and Supp(w a ) C 
Supp{wbW c ) (respectively Supp(wb) C Supp(w a w c ) and Supp{w c ) C Supp(w a Wb))- 
Identity (4.4) easily follows from cui-counting these fixed points. 

However, the general (/-case is much harder, and some negative weights remain. 
The sign of u q ((Bi, o~i); tt) comes from the cardinalities of the sets Bi and the signs 
of the permutations o^. In our proof, we successively apply five weight-preserving 
sign-reversing involutions $^ to L q {ni,n2,n 3 ), each one acting on the fixed points 
of the preceding one. $i forces 03 = Id, $2 forces S3 = 0, $3 forces o~i = Id, $4 
forces Bi = 0, and $5 forces 02 = Id, leaving B 2 arbitrary. Hence the negative 
part of (4.7) is due only to Bi- 

The final set of fixed points, Fix$5, does not contain the fixed point set (above) 
for q = 1. Instead there is a bijection from a subset of Fix&s to this set, but it 
does not preserve the powers of q. 

The five weight-preserving sign-reversing involutions ^ and their respective fixed 
points sets Fix^i are given in the next section and the complete characterization 
of F1X&5 is given by the conditions Fix.l through Fix.4, stated at the beginning of 
§6. In §6, we show that the w ? -weight of Fix& 5 is equal to the right-hand side of 
(4.5), thus establishing Theorem 3. 

5. The weight-preserving sign-reversing involutions $j. 

Let us recall that a weight-preserving sign-reversing involution (or WPSR-invo- 
lution) $ with weight function u> is an involution such that for any e ^ FixQ, 
w($(e)) = -oj(e). 

Involution $!• This WPSR-involution will kill any ((Bi, a^); tt) such that 03 is 
not the identity. 
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Remember that the cycles of a 3 are ordered by increasing minima. Find the 
greatest cycle Cj such that either this cycle is of length > 2 or it lies in the same 
block ni of tt as some other 1-cycle greater than it. If q satisfies the latter condition, 
the 1-cycle greater than a in the leftmost block tt^ of partition tt is glued to the 
end of Ci . Then, if h = ho < hi < . . . < h m = i denote the indices of the blocks 
between Wh and tti containing 1-cycles greater than Cj , these 1-cycles are moved 
from block 717^ to block TTh l _ 1 . 

For example, for ((Bi, ai); tt) G L q (9,0, 10) such that B\ = S 2 = 0, S3 = 
{10} and tt_= _( 1 )_(T,_2) (6) | (2 , 8) (5) | (_3) (4) | (5 ) (3, 9) (8 ) | (6 , 9 ) (7) (7) | (4) , we have a s = 
(l,2)(3,9)(4)_(5)(6)(7)(8) L c^=_(3 1 9), and ^((S^); tt) is given by the same 5^, 
a 3 becomes (1,2)(3, 9, 6)(4)(5)(7)(8), and ^ 

Note that the number of inversions gained in a 3 is counterbalanced by the loss 
in the statistic rs(n). Conversely, if q is of length > 2 and does not lie in the 
same block as any other greater cycles, its image is defined in the obvious way so 
that $1 is an involution. For more details, see [5]. 

Fixed points for $ x . The cycle c io is not defined if and only if a 3 contains only 
1-cycles which all lie in different blocks of tt. Therefore, 

Fix&i = {((Bi, 0"^); tt) G L q (ni, 112,713) \ 03 is the identity 
and w c contains all distinct letters}. 

Involution $ 2 - This WPSR-involution is designed to discard all ((Bi, ai); tt) G 
Fix&i such that S3 is not empty. 

Let ((Bi, ai); it) G Fix<&\ and let k = ^blocks(w). 

Denote by jo, < jo < (^3 — 1), the integer such that j + 1 = min(S 3 ). If 
S3 = 0, we let jo = 00. Likewise, denote by ji, 1 < ji < TI3, the maximum integer 
such that the 1-cycle (j x ) forms a singleton block in tt. Remember that as = Id 
and w c contains all distinct letters. By maximality, (j x ) lies in the k-th block of tt. 
Denote by j[ its contribution to the statistic rs, that is the number of (different) 
letters after the only occurrence of k in w c (and in w). If there are no such singleton 
blocks in tt, let ji = j[ = 00. 

There are two cases: jo < j[, or jo > j[. If jo < ji, &2((Bi, ai); tt) is obtained 
by inserting the 1-cycle (jo + 1) in as and by inserting the letter (k + 1) in w c at 
the (jo + l)-th position from the end of w c , leaving everything else fixed. 

For example, for ((Bi, ai); tt) defined by Bi = = B 2 , B 3 = {2,6,8} and tt = 
( 1 , 6) (5) (7) I (2) ( 1, 3, 2) (4) I (3, 5, 4) | (4) (9) | (T) | (3) | (5) , w c = 562714, j = 1, 3l = 5 and 
j[ = 2. Then the new w c in $ 2 ((Si, ai); tt) is w c = 5627184, and $ 2 ((S i7 a^; tt) is 
defined by Si = = B 2 , B 3 = {6,8} and tt = (1, 6) (5) (5) | (2) (1, 3, 2)(3)|(3, 5,4)|(4)(9)|(T)|(2)|(4)|(7). 

Note that $ 2 ((Si, ai); tt) has its j[ equal to the jo associated to ((Bi, ai); tt). 
Conversely, if j[ < jo, the image of ((Bi, ai); tt) is defined in the obvious way so 
that $ 2 is an involution. $ 2 is also weight-preserving and sign-reversing. For more 
details, see [5]. 

Fixed points for $ 2 . Fixed points correspond to the case j = j[ = 00. This 
means that S 3 = and there are no singleton blocks in tt of color 3. Therefore, 

Fix<&2 = {((Bi,ai);Tr) G Fix$i \ B 3 = and Supp(w c ) C Supp(w a Wb)}. 

Note that Supp(w c ) C Supp(w a Wb) is equivalent to the condition that the w a Wb 
is an RG- function whose maximum equals #blocks(Tr). 
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To do $3 and later $ 5 , we need to describe the contribution to the statistic rs 
of the elements of color 1 and 2 in partition tt. Let w be a word on the alphabet 
[k]. Let Wij denote the subword of w obtained by discarding letters not equal to i 
or j, 1 < i < j < k. For instance, if w = 123144124, W12 = 12112. Then we can 
write 

rs(w) = rs(wij). 

l<i<j<k 

Claim. Let w be an RG-function of maximum k and suppose w = vv' . Then v is 
an RG-function and 

rs(w) = ^2 rs ( v ij) + ^2 ls{vij) + rs{v') 

l<i<j<k, l<i<3<k, 
i$LSupp(v') iESupp(v') 

=: rs(w)\ v + rs(v'). 

Thus the contribution to the statistic rs(w) of the initial word v, rs(w)\ v , is 
indeed an interpolation between the hard statistic rs and the easy statistic Is, as was 
studied by White in [22]. He showed in particular that these specific interpolating 
statistics were (/-Stirling distributed, meaning that their generating functions over 
RG(n, k) are the (/-Stirling numbers of the second kind S q (n, k), up to a power of 
q. He provides a bijection on RG(n, k) such that the mixed statistic is sent to the 
easy statistic Is (up to a constant). More precisely, 

Lemma 4. Let S = {s\ < s 2 < ■ ■ ■ < s m } C [k]. There is a bijection '■ 
RG(n, k) — > RG(n, k) such that for any w G RG(n, k), 

m 

(5.i) Yl rs ( w ij)+ Yl ls ( w ij) = ls (^s(w))-j2( k - s j)- 

l<i<j<k, l<i<j<k, j=l 

ies ie[k]-s 



Proof. Define *j : RG(n, k) -> RG(n, k), 1 < % < k - 1 as follow: 

i) if w G RG(n, k) has a letter i to the right of the first occurrence of (i + 1), then 
the rightmost letter i is switched to (i + 1) and any (i + 1) to its right is changed 
to %. For example, *i(111212332122) = 111212332211. 

ii) if w does not have a letter i to the right of the first occurrence of (i + 1), then all 
(i + l)'s to its right are switched to i's. For example, *i(1112232) = 1112131. 
For convenience, we will set ^ : RG(n, k) — > RG(n, k) to be the identity. Now, 

given S = {s± < s 2 < ■ ■ ■ < s m } C [k], is defined as follow: 

^ s = (* fc o o . . . o * Sl ) o (* fc o o . . . o ^ S2 ) o . . . o (* fc o . . . o * Sm ). 

Note that preserves the positions of the first occurrences. For more details, the 
reader is referred to [22]. □ 

Involution $ 3 . This next involution is designed to kill any element ((Bj, Uj); 7r) 
such that o\ is not the identity. Note that since the interpolating statistics on w a 
are g-Stirling distributed, it reduces to proving the orthogonality relation 

n 

(-l) n_fc c g (n, k)S q (k, m) = 5 n , m . 

k=m 
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But this formula was deduced in Proposition 3.1 of [6] from a weight-preserving 
sign-reversing involution on appropriate pairs of 0-1 tableaux. The general idea is 
to map o\ and w a bijectively into a pair of 0-1 tableaux, using ty$ defined in the 
previous lemma and the correspondences described in §1. Then we can apply the 
WPSR-involution, essentially shifting the rightmost shortest column from one 0-1 
tableau to the other. ^s((Bi, <7j); tt) is then obtained by replacing <j\ and w a by 
the new decoded pair of 0-1 tableaux. Involution $5 will use similar ideas. 

We need only specify the bijective coding of (a±, w a ) into a pair of 0-1 tableaux. 
Let ((Bi, o~i); tt) g Fix<& 2 and let n = ni — \Bi\, k = cyc(a\) and m = max(Supp(w a )) 

i) For <7i, simply use the correspondence described in §3 to get a 0-1 tableau (f\ 
with (n — k) columns of distinct length < (n — 1). Note that inv (o"i) = inv(<pi). 

ii) For w a , we first want to reduce the interpolating statistic rs(w)\ Wa to the easy 
statistic ls(w a )- This is done by applying defined in the previous lemma to 
w a , for S = [m] \ Supp(wbW c ). We then use the correspondence described in §3 
to get a 0-1 tableau (p 2 with (k — m) columns of length < m. There is one last 
technicality: the statistic Is is sent to the non-inversion statistic on 0-1 tableaux 
(up to the constant (™)), therefore we will apply to <p 2 the symmetry involution 
exchanging non-inversions and inversions, so that for its image ip2, we have 



Note that m is not modified by the WPSR-involution applied to pairs of 0-1 
tableaux, thus insuring that the overall involution $ 3 is well-defined (the new w is 
still an RG-function) and weight-preserving. It is also sign-reversing. Details are 
left to the reader. 

Fixed points for $ 3 . At the 0-1 tableau level, the only fixed pair of 0-1 tableaux 
is (0,0), because in that case, it is impossible to move columns. But this can 
happen if and only if (n — k) = (k — m) = 0, and therefore n = k = m = n\ — \B\\i 
o\ is the identity on [m] — B±, and w a = 12 . . . (ni — \Bi\). Therefore 

Fix<&3 = {({Bi, <ji); tt) G Fix<& 2 I ci is the identity and w a = 12 . . . (ni — |i?i|)}. 

Involution $4. This involution is the simplest. Its task is to eliminate elements 
((B i: <Ti); tt) such that B\ 7^ 0. 

Let ((Bi, o~i); tt) G Fix^s and let io be the smallest integer, 1 < zo < ni, such 
that either i G Bi, or the 1-cycle (i ) forms a singleton block in tt. Then if i G Bi, 
insert it as a 1-cycle in o\ and as a singleton block in tt, and vice-versa. 

For example, if B x = {2},S 2 = B 3 = 0, and tt = (1)(T, 3)|(3)|(1, 2) (2), then 
io = 2 and the image of ((Bi, cr^); tt) under $ 4 is Bi = 0, B 2 = S3 = 0, and 
tt = (1) (I,3)|(2)| (3) I (1,2) (2). Details are left to the reader. 

Fixed points for $4. 

Fix<&4 = {((Bi,<ji);TT) G Fix^s \ B\ = and Supp(w a ) = [ni] C Supp(wbW c )}. 

Involution $ 5 . This final WPSR-involution will annihilate the remaining ((Bi, cr^); tt) 
such that a 2 is not the identity. It is the only one using the hypothesis 123 > n\ > n 2 . 
The principle of the involution is similar to $3: we will reduce the problem to find- 
ing an involution for the easy statistic Is. 
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Let ((Bi, (Tj); n) G Fix<&4, and let ^blocks(ir) = ns + s. First, encode 02 as a 0-1 
tableau </? with {rii — \Bi\ — cyc(o~2)) columns of distinct lengths < (77-2 — I-B2I — 1) 5 
using the correspondence described in §3. Note that inv(a2) = inv(ip) and that the 
shortest column of <p> is of length at most cyc(o2)- 

For Wb, we reduce the interpolating statistic rs(w)\ WaWb to the easy statistic Is by 
applying defined in Lemma 4 to w a Wb, with £ = [773 + s] — Supp(w c ). Note that 
since w a = 12 . . .ri\ and preserves first occurrences, ^ s(w a Wb) = w a Wb for some 
word Wb — b\b<2. . . . b&- Note also that we must have {ni + 1, ... , n 3 + s} C Supp(wb) 
(because w a Wb has maximum (77.3 + s)). 

For example, if ({Bi, c^); n) G Fix<&4 is defined by B\ = B^ = B% = 0, and 
^ = ( 1 ) (2) I (2) (3) (5) j (3) (2) I (4) (4) j (5) (5) (T) | ( 1 ) (4) (3) , we have w a = 12345, w b = 
61265, w c = 53642, and 02 = (1)(2)(3)(4)(5). Then cr 2 corresponds to the empty 
0-1 tableau (p = 0, and we successively compute S = [6] — Supp(53642) = {1}, 
y {1} (w a w b ) = 1234566154, and w b = 66154. 

Let io denote the length of the shortest column in ip, 1 < io < cyc(p2)- If </? = 0, 
let io = 00. Likewise, let ho denote the smallest integer, 1 < ho < cyc(<j2), such 
that bh < ho- If no such bi exists, set ho = 00. 

There are two cases: io > ho or io < ho- If io > ^0, then delete bh from the 
word Wb and add a column of length (ho — 1) to <p, with a 1 in position bh , from 
bottom to top, thus obtaining a new pair (w' h , tp'). Since the letter removed from 
Wb is at most equal to (cyc(a2) — 1) < (112 — I-B2I) < (^1 + 1), w a w' b is still an RG- 
function of maximum (ns + s), and the new io associated to ip' is equal to (ho — 1). 
$5 ((Bi, (7j); 7r) is then obtained by applying to w a w' b and by decoding the 0-1 
tableau ip'. 

In the above example, io = 00 and ho = 3. Hence ip' = 



1 



(corresponding to 



the new permutation 02 = (1)(2, 3)(4)(5)) and w' b = 6654. From ^ ^(w a Wb) = 

123456165, we get $ 5 ((^, O; tt) equals B x = B 2 = B 3 = 0, and tt = (1) (2, 3) j (2) (5) | (3) (2) j (4) (4) j (5) (5) ( 

(1)(4)(3). 

If io < ho, the image of ((Bi, a^); it) is defined in the obvious way so that $5 is 
an involution. The proof that $5 is weight-preserving and sign-reversing is quite 
straight-forward, and the details will be left to the reader. It remains to show 
that $5 is well-defined. Remember that if ((Bi, cr^); tt) G Fix<&4, we must have 
Supp(w a ) C Supp(wbW c ). We have to show that $5 preserves this property. What 
complicates matters is the application of and ^g 1 to the RG-functions w a Wb- 
In Lemma 5, we explicitly find the set of images w a Wb (which we will denote by 
W(S)) of all possible w a Wb under ^5. We will then show that the deletion or 
insertion of a letter whose value is strictly less than its position in u>b yields new 
RG-functions w a w' b which remain in the set W(S). 

Fix ns > n\ > 112 > 0, < t < 712, and < s < n\ + 712 — ns. Let S C [n^ + s] 
such that \S\ < s, and fix w a = 12 . . .n\. We denote by 



W(S) = {wb I w a Wb G RG(ni + n 2 — t,n 3 + s), and 
[ni] C ([ti 3 + s] - S) U Supp(w h )}, 

W(S) = {w b I ^ s (w a Wb) = w a w b for w b G W(S)}, 
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and 



w a W(S) = {w a w b | w b G W(S)}, 
w a W(S) = {w a w b | w b G W(S)}. 

In particular, when |,S| = s, if w c is a word containing the letters in ([n3+s] — S) in 
any order, with no repetition, and (B 2 , a 2 ) is a partial permutation of {2} x [n 2 ] with 
cyc(a 2 ) = n 2 — t, W(S) contains all possible words Wf, such that w = 12 . . .n±WbW c 
is the RG-function associated to some ((-Bi, o~i); n) G Fix<&4 having these fixed 
(£? 2 , cr 2 ) and iy c . 

Lemma 5. (characterization of W(S)) Let S C [77,3 + s] sitc/i IS"! < s. The 
set W(S) depends only upon the cardinality j = \S fl [ni]|. More precisely, we have 

(i) W(S) = W(Sn[n 1 }), 

(ii) If j = 0, W(0) = W(0), and Wb G W(0) has the following form: 
(5.2) 

Wb = (ni + 1) (ni + 2) . . . (713 + s — 1) *_^5 ( n 3 + s ) • 

entries <n 1 < (ni + 1) < (n 3 +s — 1) < (n 3 +s) 

(mj if j = 1, t/ien VF({z}) = is obtained from W(0) by keeping only the 

words Wb of the form (5.2) such that one of the stars * is set to its maximum 
and the maximum value of all the stars to its right is lowered by 1. So any Wb 
has the form 

Wb = ( n i + 1) ( ni + 2) . . . (ni + h) ( n i + h) 

entries < ni < (rti + 1) < (n 1 +h) <(tii+/i- 1) 

(5.3) 

(ni + h + 1) . . . (ri3 + s — 1) ,*_^5 ( n 3 + s ) ij^i3 • 

<(n 3 +s-2) <(n 3 +s-l) 

(wj #j > 2, i/ien = W({1,2, ... , j}) is obtained from W '({1,2, .. . ,j-l}) 6?/ 

£/ie same construction as the one described in (Hi). 

Proof. 

(i) . First we show that W( S) = W(Sf][n 1 }). From the definition of W(S), it is clear 
that W(S) = W(S n [ni]). Moreover, if S = {s 1 < . . . < Sj < s j+1 < . . . < s n }, 
where Sj < n\ and Sj+i > n\, since \I/s\[ ni ] is a bijection on RG(ni + n 2 —t,ri3 + s), 
preserving first occurrences and leaving all letters < ri\ fixed, we must have 

y S \[ ni ](waW(S)) = w a W(S). 

Therefore, 

w a W(S) = * s (w a W(S)) = *sn[ ni ] ° *s\[ ni ](waW(S)) 
= * SnW (^(5n[m])) = w a W(Sn H). 

(ii) . If j = 0, \l/0 is the identity map and 

W(0) = W(0) = {w b \12... niw b G RG(n! + n 2 -t,n 3 + s)}, 
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in which typical elements (tails of RG-functions) are given by (5.2). 
(iii). If j = 1, suppose S = {z}, 1 < i < n\. Then 



W(S) = {wb | w a Wb G RG(ni + n 2 — t, ns + s), and i G Supp(wb)}- 

Let Wb G W(S) and suppose the rightmost occurrence of i lies in position p of Wb, 
between the first occurrence of (ni + h) and the first occurrence of (ni + h + 1). 
Thus w a Wb has the form 

w a Wb =12 . . . ni ^_^3 ( n i + 1) ( ni + 2) • • • (ni + h) ^-L^3 

entries < ni < (ni + 1) < (rti +h) position < (rti+/i), 

(ni+p) entries j^i 




<(n 3 +s-l), < (n 3 +s), 

jti & 

Apply ^{i} = ^n 3 +s ° ^n 3 +s-i • • • ° to w a Wb- The last occurrence of z in 
w a Wb (in position (ni + p)) lies to the right of the first occurrence of (i + 1) (case 
(i) in the definition of \I/ m ), so it is changed to (i + 1) by and any (i + 1) to its 
right is changed to i. Thus the last occurrence of (i + 1) in ^i(w a Wb) now appears 
in position (ni + p), again to the right of the first occurrence of (i + 2). So all 
(i + 2)'s to its right are changed to (z + l)'s by the (i + 1) in position (ni +p) 

is switched to (i + 2), and every other letter remain fixed. 

The same argument applies until we reach \E f ni+ / l . At this point in \l/ ni+ / l _i o 
. . . o \l/j(t(; a t(; b ), there is a (ni + h) in position (ni +p) and no occurrence of (ni + h) 
to its right. This means that there are no letters (ni + h) to the right of the first 
occurrence of (ni + h + 1) (case (ii) in the definition of ^ m ). Hence ty ni+ h changes 
every occurrence of (ni + h + 1), except for the first one, to (ni + h)'s, and fixes 
everything else. Once again in the RG-function obtained, there are no occurrences 
of (ni + h + 1) to the right of the first occurrence of (ni + h + 2). It is clear that 
by applying successively * ni+ ^ + i, . . . , ^ns+s respectively, we will get ^{i}(w a w b ) 
exactly of the form (5.3). This shows that the set defined in (iii) is equal to ^({i}). 
Note that the definition of the set ^({i}) is independent of the actual value of z, 
so W({i}) = W({1}). 

(iv). The proof is an easy induction based on the proof of (iii). Note that if 
>S = {si < S2 < ■ ■ ■ < Sj},Sj < ni, the positions of the last occurrences of 
si, S2, ■ ■ ■ , Sj respectively in w a Wb correspond exactly to the positions of the stars 
successively fixed to their maximum in ty s (w a Wb). □ 

We can show now that $ 5 is well-defined. 

Let Wb G W^({1, 2, . . . , j}). The letters of Wb can be divided into two categories: 
the fixed letters (first occurrences of (ni + 1) up to (ns + s), and j stars that were 
fixed to their maximum in the construction described in the preceding lemma), and 
the free letters (corresponding to stars in the description of Wb in Lemma 5). So in 
order to be in 2, . . . , j}), a word Wb must have (713 + s — rt\ + j) fixed letters 

(appearing in some fixed relative order), and possibly some free letters, depending 
on its length. 

On one hand, note that the fixed letters of w>b are always greater or equal to 
their positions in w b - Indeed, we have already seen that the first occurrences were 
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necessarily greater than their position p ((77.3 + s) > {n\ + 1) > n 2 > p). As for the 
j stars fixed to their maximum, the way to minimize their value in the construction 
of Lemma 5 is to fix them successively by increasing order of their positions. Then, 
if they all lie before the first occurrences of (ni + 1) up to (77-3 + s), the j-th star 
fixed will have minimum value (77 1 — j + 1), and the rightmost position where it 
can be located is, for example, the one in the following word: 

Wb = **...* n\ (m — 1) . . . (ni — j + l)(ni + l)(m + 2) . . . (773 + s). 

entries< ni 

But from the relations 77,3 > 77-2, t > 0, and j < s, we deduce that its position p, 

P = \w b \ - (n 3 + S - Tlx) = 77i + (n 2 - 77 3 ) - t - S < Til ~ j + 1- 

Therefore in that case, all fixed stars are greater or equal to their positions. More 
generally, if a fixed star is rather located to the right of a first occurrence, its value 
is increased by one, so the letter remains greater or equal to its position. 

On the other hand, note that the allowed maxima for the free letters are also 
greater or equal to their positions in w^. The same type of argument (with same 
inequalities) applies. Details are left to the reader. 

Now, the "involutive step" of $5 was to add or to delete a letter from Wb, and 
this letter had the property of being strictly smaller than its position in Wb- 

If the involutive step deleted a letter from u>b (w a w b G RG(ni + 712 — t, 77.3 + s)), 
then it had to be one of its free letters because the fixed ones are greater or equal to 
their positions. Therefore the new w' b obtained is in the set 2, . . . , j}) (with 

w a w' b G RG{n\ +77-2 — (t + 1), 773 + s)). Likewise, if the involutive step added a 
letter to Wb, the new letter is in the right range to be considered a free letter, and 
the fixed letters (and their relative order) are not modified, so the new w' b is in the 
set 2, . . . , j}) as well (with w a w' b G RG{n\ + 77,2 — (t — 1), n 3 + s)). 

Fixed points for $ 5 . The fixed points of $5 correspond to the case i = h = 00. 
Clearly, 

Fix<& 5 = {((Bi, Oi);n) G Fix$ 4 \a 2 = Id and for S = [#blocks(ir)] — Supp(w c ), 
the word Wb in ^s( w a w b) = w a Wb has its z-th letter > z,Vi}. 

6. Combinatorial evaluation of L q (C ni (x)C n2 (x)C n3 (x)). 

An expression of L q (C ni (x)C n2 (x)C n3 (x)) can now be computed by w g -counting 
of the remaining fixed points Fix& 5 . More precisely, ((Bj, Uj); 7r) G Fix& 5 if and 
only if 

Fix.l B 1 = B 3 = 0, 

Fix.2 <n = Id for i = 1,2,3, 

Fix. 3 w a (respectively w c ) has all distinct letters and Supp(w a ) C Supp(wbW c ) 
(respec- 
tively Supp(w c ) C Supp(w a Wb)), 
Fix. 4 for S = [#blocks(ir)]—Supp(w c ), the word w b = bib 2 ■ ■ ■ 6 n2 _|s 2 | in ^s(w a Wb) = 

w a Wb has all b^s > i, where was defined in Lemma 4. 
Clearly, for such elements, the weight (as was defined in (4.8)) reduces to 

(6.1) W q ((B h (Ji) : 7r) = (_l)|S 2 | 5 int»(B 2 )+r a (7r) a |B 2 H-#Wocfc fl W_ 
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By cu 9 -counting this fixed point set, we will show that 



(6.2) 



ni+n 2 —n 3 I s 

L q (C ni (x)C n2 (x)C n3 (x)) = £ EE^Hl'^W^ 

1=0 s=0 j=0 



n 2 
l-s 



m 




n 3 - ni + s 




n 2 - / + s 


_ 3 _ 


q 




q 


n 3 - ni + s 



bV- q [ni - j]\ q 
[n 3 -n 2 + l]l t _ 



where 



s-J 
2 



ni + n 2 - n 3 - / 
3 



- (s-j)(n 3 - ni +j) 



fn 2 — I + s\ 

+ 1 2 I + (n 3 - n x + s){n 3 - n 2 + I) + j{n 3 - n 2 + I). 

Evaluating the s-sum by the g-binomial theorem (which has a simple bijective 
proof) gives the right-hand side of (4.5) and thus Theorem 3. 

The main difficulty here is to transpose the condition Fix.4 into the Wg-counting. 
Using Lemmas 4 and 5, we will see that this corresponds to the (/-counting of some 
special sets of RG-functions according to the statistic Is, which is the object of 
Lemma 6. 

Let us first group the elements of FixQ*, by powers of a. The power of a ranges 
from a minimum of n 3 (expressing the fact that w c has n 3 distinct letters) to a 
maximum of (ni + n 2 ) (being the maximum value of {m8ix(Supp(w a wi ) )) + l-E^I})- 
Now, 



(6.3) 



ni+n 2 —n 3 



L q (C ni C n2 C n3 ) = £ a"3+^(_iy-V ! r) 



1=0 



s=0 



n 2 
l-s 



q ir 



where the last sum ranges over all partitions tv corresponding to ((Bi, o^); n) G 
Fix&s such that ^blocks(rr) = (n 3 + s) and B 2 is any fixed subset of {2} x [n 2 ] of 
cardinality (I — s). The s-sum is the generating function for the subsets B 2 , as was 
established in §2. But 



(6.4) 



rs(Tv) = rs(w) = rs(w c ) + rs(w)\ WaWb 

= rs(w c ) + ls(w a w b ) - ^ 



(n 3 + s - u). 



uE([n 3 +s] — Supp(w c )) 



Note that for any fixed set Supp(w c ), there are no constraints on the positions of 
the letters in w c , so rs(w c ) is simply the number of inversions of the word w c , whose 
distribution is mahonian (i.e. the generating function equals [n 3 ]!q). From Lemma 
5, we also know that the possible choices for Wb only depend on the cardinality j of 
the set ([n 3 + s] — Supp(w c )) fl [ni], not on the actual set Supp(w c ) itself. Hence, 
if we let 



FixW(j) = {w b | w b = b 1 ...b n2 -i +s E W({1,2,... ,j}) and h > i,Vi}, 
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where W({1,2, . . . was characterized in Lemma 5, we get that the last sum on 
the right-hand side of (6.3) equals 



n J=° w b £FixW(j) 



ls(w a w b ) 



(6.5) 



gj'(-«3-« + l)+(5) 



Tlx 

3 



-( S 2 J )-( s -j')("3-ni+j) 



n 3 — ni + s 
s - 3 



Finally, we show that 

Lemma 6. Ifw a = 12 . . . n\, Wb = b\ . . . b n2 -i +s and max (Supp(w a Wb)) = n 3 + s, 
then 
(6.6) 



w b EFixW(j) 

where 



qis(w a w b ) _ yA 



n 2 -l + s 
n 3 - ni + s 



[3?. q [ni - J]l q 
[n 3 - n 2 + /]! ( 



n\ + n 2 - n 3 - I 
3 




Proof. Note that the statistic Is of any RG- function is just the sum of the values 
of the letters minus one, so 



(n 2 — l+s 
7 i=i 



ls(w a w b ) 



To visualize more easily where the various factors of (6.6) come from, let us 
encode as a 0-1 tableau <p in the following manner: start with a (n 3 + s) x (n 2 — 
/ + s) rectangular Ferrers diagram. Fill it with a 1 in position j (from bottom to 
top) of column i if hi = j, and with O's elsewhere. For example, if (n 3 + s) = 8, 
ni = 6 = (n 2 — l + s) and Wb = 175787, is the 0-1 tableau on the left of Figure 3. 
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Figure 3: Encoding of Wb as a 0-1 tableau 

Obviously, we have XX^* ~ -0 = inv(ip). Note also that the O's in the shaded 
staircase shape of f in Figure 3 always count as inversions, expressing the fact 

~ ( ri2 — I ~\~ s \ 

that hi > i. They account for the factor q\ 2 ' in (6.6). We can therefore drop 
them from cp without loss of generality, and compute the inversion number of the 
reduced 0-1 tableau dp. We now use Lemma 5 to characterize the possible fillings 
of (p according to j. 
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Case 1: j = 0. From Lemma 5 (ii), Wf, G W(0) simply means that it is the tail 
of an RG-function. Thus the only restrictions on are that the first occurrences 
of (ni + 1), (toi + 2), . . . , (n 3 + s) appear in the right order. 

If we set x = (tt-3 + s), y = ni, and z = (p,2 — I + s), in the context of 0-1 
tableaux, we want to g-count all 0-1 tableaux Cp with z columns of lengths x, (x — 
1), . . . , (x — z + 1) respectively, such that when we look at the top (x — y) rows of 
<p from left to right, the leftmost 1 in any row must always occur before the ones in 
the rows above it. Grouping the tableaux according to these leftmost occurrences 
of l's, we get "typical" 0-1 tableaux 0typ, corresponding exactly to the typical 
words Wb described in (5.2) of Lemma 5. For instance, the typical 0-1 tableau (pt yp 
containing our previous example is illustrated in Figure 3 (stars * correspond to 
possible positions of l's). Carrying out the (/-counting, observe that 

1.1 Each column containing a number m of stars contributes a factor [m] q to the 
(/-counting of inversions. No matter which (x — y) columns are chosen to be 
first occurrences of upper l's, the number of stars in the remaining columns is 
y, (y — 1), . . . and (x — z + 1) respectively, contributing to an overall factor of 

[y% [ni]\ q 
[x-z]\ q [n 3 - n 2 + l]\ q 

1.2 The 0's below the leftmost occurrences of upper l's (shaded in Figure 3) form a 
partition \i with (x — y) parts of length at least (x — z) and at most y, determined 
by the positions of the first occurrences. Summing over all possible choices, it 
contributes a factor 

n 2 - I + s 



(x-y)(x-z) 



Z 

x-y 



nz — n\ + s 



J q 



q 

Case 2: j > 1. Recall that Lemma 5 (iii) and (iv) provides a method to construct 
all the elements of 2, . . . , j}) uniquely from W(0). In the 0-1 tableau con- 

text, if we extract only the cells filled with stars in (p typ (hence obtaining a tableau 
iptyp with (m + 77-2 — ns — I) columns of lengths ni, (ni — 1), . . . , (77-3 — 77-2 + I + 1) 
respectively), the manipulation described in Lemma 5 (iii) corresponds to replacing 
the top star of a column by a 1, and all the top stars to its right and the stars below 
it by a 0. Repeating this procedure j times and reinserting the columns of ifj t y P in 
(fityp yields to "typical" 0-1 tableaux that correspond to the elements of FixW(j). 
For example, Figure 4 shows the above manipulations on the third and the first 
columns respectively of the stars extracted from (p tyv of Figure 3. 
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Figure 4: Manipulations of Lemma 5 in the context of 0-1 tableaux 
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Proceeding to (/-counting, the part (1.2) of case 1 is left unchanged and the part 
(1.1) is replaced by the contribution of the different choices of iptyp- But observe 
that 

2.1 All the 0-1 tableaux ifjt yp such that a star has been changed to a 1 in columns 
ci, C2, • • • and Cj contribute to [j]l q times the g-counting of the 0-1 tableaux ifj typ 
such that this procedure was done in increasing order of the Cj's. Therefore, we 
can restrict to this latter case. This explains the factor [j]\ q in (6.6). 

2.2 It is not hard to see that in that case, we are g-counting all 0-1 tableaux if) 
containing (m + ri2 — n$ — I) columns of lengths n\, {n\ — 1), . . . , (77-3 — 77,2 + 1 + 1) 
respectively, such that when we look at the top j rows from right to left, the 
rightmost 1 in any row has to occur before the ones in the rows above it. There 
is a simple weight-preserving bijection between this class of 0-1 tableaux and the 
one that was g-counted in case 1, for x = n\, y = (n\—j) and z = (771+77,2—77,3—/) 
(this class is defined by interchanging "left" and "right"). Given ip in the first 
class of 0-1 tableaux, just leave all the l's below the j-th row fixed and "reverse 
the order" of the l's in the top j rows, within the columns where they appear. 
Figure 5 gives an example of this for j = 2. 
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Figure 5: Bijection between two classes of 0-1 tableaux 

This is clearly an involution that preserves the number of 0's below l's. There- 
fore, we can simply use case 1 to compute the g-contribution of the V>'s. We obtain 



J(n 3 -n 2 +l) . 



[m -3Y-c 



[n 3 -n 2 + /]!< 



n\ + ri2 — 773 — I 
3 



□ 



Finally, putting together Lemma 6, identities (6.5) and (6.3) yields identity (6.2), 
thus completing the proof of Theorem 3. 

Note that if we take 77,2 = and apply $1 and $2 to L q (ni, 0, 77-3) (assuming 
^3 > 77.1)5 the set Fix<&2 is easily seen to be empty unless n\ =773, in which case it 
can be proven to be Wq-counted by the right-hand side of (4.2) (n = ni,m = 773), 
thus proving orthogonality and Theorem 2. These results can also be obtained by 
applying directly $3 and $4 to the set L q (ni, 0, 773), assuming this time tt-i > 773. 
We also have a weight-preserving sign-reversing involution proving orthogonality 
when the g-statistic for the moments is taken to be lb instead of rs, but we do not 
know how to generalize it to the linearization problem. 

Corollary 7. Let n\ > n 2 > . . . > 77/-. The coefficient of the lowest power of a, 
a ni in L q (C ni C n2 . . . C nk ) is a polynomial in q with positive coefficients. 
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Proof. The proof of Theorem 3 can be generalized to a product of k g-Charlier 
polynomials, any additional color being treated as was color 2, the middle color. It 
is easy to see then that the fixed points contributing to the lowest power of a must 
have all Bi — 0, and therefore have all positive weights. □ 

Corollary 8. Let n 3 > n\ > n 2 . The coefficient of a ni+n2 ~ l in L q (C ni C n2 C n3 ) is 
equal to (q — i) n i+ n 2-n 3 -2i n mes th e coefficient of a n3+l , for < i < [(ni + n 2 — 

Our proof of Corollary 8 is analytical, but we would like to have a combinatorial 
explanation of this "symmetry" property. 

Note that Fix&$ is not an optimal set of fixed points, in the sense that there are 
still some terms that cancel each other when we proceed to Wg-counting of FixQ$. 
For example, for n\ = n 2 = ns = 2, the two elements of Fix<&$ such that B 2 = {2}, 
w = 12121 and B 2 = 0, w = 123123 have weight — a s q 3 and a 3 q 3 respectively. 
However, we do not believe that an attempt to reduce Fix<&$ would be worthwhile. 

Corollary 9. Let n x > n 2 > . . . > n k . If q = 1 + r, L q (C ni C n2 . . . C n J is a 
polynomial in r with positive coefficients. 

7. The classical g-Charlier polynomials. 

We contrast the results of the previous sections with those for the classical q- 
Charlier polynomials [11, p. 187] 

(7.1) c n {x-a-q)= 2 0i(<T n , x; 0; q, -q n+1 /a). 
The monic form of these polynomials, cc n (x;a;q) satisfies 

cc n+1 (x; a; q) = (x - b n )cc n (x; a; q) - A n cc n _i(>; a; q), 

where 

A n = -a^- 2n (l - q- n )(l + aq~ n ), b n = aq- 1 ^ + q~ n + aq~ 2n - aq~ n . 
A calculation (see [11, p. 187]) shows that the moments for these polynomials are 

n 

= Y[(l + aq~ l ). 

i=l 

We need to rescale x and a so that b n and A n are (/-analogues of a + n and an 
respectively. If we put x = 1 + z(l — q), and multiply a by (1 — q), and call the 
resulting monic polynomials C n (z; a; q), the explicit formula from (7.1) is 

n 

(7.2) C n (z;a;q) = q- n2 J2 

fc=0 

The three term recurrence relation coefficients are 

(7.3) b n = q- n [n) q (l + a(l-q)q- n ) + aq- 1 - 2n , X n = aq 1 - 3n [n] q (l + a(l-q)q- n )- 
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{-a) n - k q( h t 1 ) 1 [[(q i z-[{\ q ) 



A calculation using the measure in [11, p. 187] gives 

n 

(7.4) v^J^Q-^SihinJW- 

i=i 

Again we find (/-Stirling numbers for the moments. Zeng [24] has also derived (7.2) 
and (7.3) from the continued fraction for the moment generating function. 

We see that the individual terms in (7.3) do not have constant sign. This means 
that the Viennot theory must involve a sign-reversing involution for its combinato- 
rial versions of (7.3) and (7.4). Nonetheless we can give combinatorial interpreta- 
tions of (7.2) and (7.4), but have no perfect analog of Theorem 3. 

Acknowledgement. Theorems 1 and 2 were originally found in joint work with 
Mourad Ismail. 
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